Localization transition of random copolymers at interfaces 
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We consider adsorption of random copolymer chains onto an interface within the model of 
Garel et al. Europhysics Letters 8, 9 (1989). By using the replica method the adsorption of the 
copolymer at the interface is mapped onto the problem of finding the ground state of a quantum 
mechanical Hamiltonian. To study this ground state we introduce a novel variational principle for the 
Green's function, which generalizes the well-known Rayleigh-Ritz method of Quantum Mechanics 
to nonstationary states. Minimization with an appropriate trial Green's function enables us to find 
the phase diagram for the localization-dclocalization transition for an ideal random copolymer at 
the interface. 



PACS numbers: 61.41. +e, 05.40. +j, 03.65.-w 

The presence of copolymers at interfaces between two immiscible fluids is crucial to such processes as emulsion 
stabilization, wetting, microemulsion formation and reinforcement of polymer-polymer interfaces |^- ||l0[| . For oil- 
water interfaces, the polymers that are used in these applications are amphiphilic in nature. While one component is 
soluble in the oil phase, the other component is soluble in water. The difference in solubilities drives the copolymers 
to adsorb at the interface between the two phases. The localized copolymers can stabilize the interface in the sense 
that they significantly reduce the surface tension. The study of random copolymers has been motivated in recent 
years by relevance of these materials for both biological and technological applications. Moreover, the properties of 
these simple random systems may be important in understanding of much more complex systems such as proteins 

It is well-known that adsorption of long polymer chains on surfaces or interfaces is related to the bound state 
problem of a Quantum Mechanical (QM) particle in a potential well The Green's function G(r,r';s,s') of a 

chain in d spatial dimensions with the monomer s at position r and the monomer s' at r' is a fundamental quantity in 
the statistical mechanics of polymers. It obeys in the presence of an external potential ^(r, s) the following equation 

with the condition: G(r, r', 0) = d{r — r') and I being the statistical segment length of the polymer. Eq.(|l|) is related 
to the Schrodinger equation by using the replacement: s it, P /{dkT) kT fi [Q. Considering an 

external potential independent of the monomer species along the chain and given by V{x)/kT = —uS{x) in one spatial 
dimension, a localized part of the solution exist and can be written as Gioc{x,x' ,N) — fc exp(/^A^fc^/2)) exp(— fc|x| — 
fc|a;'|). It becomes the only relevant contribution for TV = |s — s'| ^ oo. Here the localization length ^ = 1/fc is given 

by e = l/u. 

However, random copolymers utilized in order to reinforce polymer-polymer interfaces do not adsorb according 
to the above scenario. The external potential in that case has an arc-length dependence, which in QM picture 
corresponds to a time dependent potential, so that an appearance of a bound state is not so obvious. In contrast to 
the example considered above, the individual monomers belonging to a random copolymer are not a priori attracted 
by the interface. The localization effect appears to be due to the feed-back of the interaction of larger parts of the 
chain with the interface. Garel et al. [|l6| introduced a simple model describing the behavior of an ideal copolymer 
chain in the presence of an interface between two phases by the following Edwards functional integral 

V[x{s)]exp{-— ds{dx/ds f+w d<(s)sgn(x(s))}, (2) 

where C(s) is a random Gaussian variable describing the heterogeneity of the copolymer. Due to the fact that there 
are no interactions between the monomers, the coordinates y{s) and z{s) along the interface separate, so that the 
problem becomes one-dimensional. In the model of Garel et al. ]l6[ | the distribution function of the random variable 
C(s) reads for the discrete version of the polymer chain (x(s) — s- Xi, ^ di i = 1, ■■■,N with N being the degree 
of polymerization of the chain) P(Ci) = (27rAo)~^/^ exp(— ((^ — (o)^/2Ao), where Co describes the asymmetry in the 
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composition of the copolymer. The average over Q within the rephca method ||T^ leads to the following Replica- 
Hamiltonian 

n n ^ n 

Hn ^ -D^\/l~ wCo^sgn{xa) ~ -^Aqw'^ X! '5gn(xa)sgn(x/3), (3) 

a— 1 a— 1 Q,/3— 1 

where n is the number of replicas and D is given by D = So far, besides Imry-Ma type arguments |^,|| there is 

no analytical treatment enabling one to find the ground state of Hamiltonian 7J„ , which is related to the localization of 
the random copolymer. The Hartree method applied in Ref. fails to describe the localized state. The problem at 
hand it also of general interest since it is a further example where frozen-in disorder (here chemical disorder) induces 
localization of the chain. 

In this Letter we present a novel variational principle in order to study the ground state of the Hamiltonian 
given in Eq.(^). In contrast to the well-known stationary Rayleigh-Ritz (RR) method |2^, where the function of 
interest is the wave function of the problem, we are interested in a variational principle for the full Green's function 
(GF). While the GF is a dynamic quantity (in the QM picture), the latter is a generalization of RR extremal principle 
to nonstationary states. We note that our variational principle for the Green's function is not restricted to treat the 
Hamiltonian iJ„, Eq.(^, but it is in fact a general method of quantum mechanics and in particular is a promising 
method for treating polymer problems. 

In order to introduce the variational principle we start with the definition of the GF written in the form 

- G-\z) + G^\z) + = 0, (4) 

where Go is the Greens function without interaction, Gq^{z) — z + Hq, Hq is the unperturbed part of the Hamilto- 
nian, z is Laplace conjugate to the chain's length (time in QM), and Hi is the interaction part of the Hamiltonian. 
Considering Eq.(||) as a stationarity condition, 5T{G)/5G = 0, of a functional T{G), we obtain 

T{G) = -trln(G) + trGf^^G + trH.G. (5) 

Eq.(^ defines an extremum functional, the stationarity condition of which yields the exact Green's function. We now 
present a refined version of this principle. Iterating Eq.(^) we arrive at 

^ (1 - H,Go)G-' + Go 1 - H,GoH, = 0. (6) 

Considering the latter as a stationarity condition for a functional T'{G), we get 

T'{G) = -tr(l - H,Gq) ln(G) + trG^^G - tiH.GoH.G. (7) 

For practical purpose it may be useful to start with T' instead of JF. We remark that T{G) and especially T'{G) are 
similar to the generating functional of the 2nd Legendre transform in statistical physics and in field theory [0- [^ . 
It is amazing that the method based on (|5[ 0) has been discovered only now. 

We now will apply the above method to find the ground state of an asymmetric localization potential given by 
U{x) = —u5{x) -f x(^{^)^ where 6{x) is the step function. Assuming ground state dominance (GSD) iQ we choose 
the trial Green's hmction as 

G(ki,k2;x,x' ,t) — -— !— ^ exp(— £fct)(exp(— /cia;)f?(a;) -|- exp(fc2x)0(— x)) 
fci -I- k2 

{ex-p{-kix')e{x') + exp{k2x')e{-x')), (8) 
where Ek = — and fci and k2 are fit parameters. Inserting G into T{G) yields 
T(i, h \ w , ^ , + Dkik2 u 2kik2 x 

T{ki,k2) = \n{z + ek) + — , +— — , , , ■ 9) 

z + £k z + Ek ki + k2 z + Ek ki + k2 

Solving equations associated with stationarity conditions 

dT/dki = dJ'/dk2 = (10) 

we get the fit parameters ki and k2 as 
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^1^0 n+~' ^2 = 7^7^--. 11 
2, D u Z D u 

The latter coincide with those obtained from the exact solution of the Schrodinger equation in the asymmetric 
potential. For disappearing asymmetry (x — 0) the ground state will become symmetric (fci = k2)- The delocalization 
transition is defined by the condition fc2 = 0. We note that the minimization procedure based on (|^) gives the same 
result as the Rayleigh-Ritz approach. We can also consider Sk in Eq.(|9|) as a free parameter and minimize J-{£, fci, 
given by Eq.(^ with respect to e, ki and ^2- The stationarity condition with respect to e fixes e as the expectation 
value of the Hamiltonian associated with (^. The stationarity conditions with respect to ki and ^2 coincide exactly 
with those obtained by using the Rayleigh-Ritz method. 

We now turn to the adsorption of a random copolymer chain, or what is equivalent to the study of the ground 
state of the Hamiltonian 7i„ , Eq. (^) . In the framework of the variational principle based on the functional J^, which 
is given by Eq.(||), the simplest n-replica trial Green's function is 

2k k " 
G(fci, fc2; X, x' , t) = (- — —tT exp{-nekt) TT {{exp{-kiXa)d{xa) + exp{k2Xa)d{-Xa)) 
ki+k2 ^ 

a—l 

{eM-kix'jeix'J + eMk2x'J9{-x'J)) (12) 

with Ek = ^Dk2 in the case Co > 0. The energy Sk is negative due to the same reason as in the case of the localization 
of a QM particle. Notice that the trial function in Eq.(ll^) is a product over the one replica Green's functions. This 
is due to the circumstance that the main effect of the randomness here is expected to result in an attraction of the 
copolymer to the interface. In the case when the randomness would result in self-interactions between the monomers, 
a product over replica pairs would be more appropriate. Eq.([l^) can be considered as the first term in the expansion 
of the n-replica Green's function over one-replica, two-replica, ... Green's functions. The evaluation of J- gives 

-nek + nDkik2 wCqu k2 - ki n ^ 1 n(n - 1) ^ 1 M-ki ^ 
T{ki, k2) ^ ln(z + nsk) + - -Aqw — Aqw — ( ) . 

z + nek z + nek ki + K2 z z + nek 2 z + nek ki + k2 

(13) 

After some algebra the stationarity conditions ( |lO| ) associated with (13) simplify in the limit n = to 



k-^ + 3 k^ki + 3 fc2 fcf + fcf + 2 C fc2 + 2 C fci - 2 A fc2 + 2 A fci = 0, (14) 

kj + 3 klki - fc^C + 3 klkl - 2 A fci fc2 + fc2 fci C fc? = 0, (15) 

where the quantities A = Aqw'^ /D and C = (^qw/D are introduced. The solution for ki and k2 is obtained as /c2 — aki, 
where a is given by a = (C + 2 A + 5)/(2 A-3 C) with S = a/4 + 2 A ( + 4 A^. The quantities fci and ^2 are obtained 
from Eqs.(||-^ for C > as 



(-2C + 4A + S')^/^ 



(16) 



.2^v/2(C-.2A-^.)(^^l^i^±M^ (17) 
(-2C + 4A + 5)'/' 

The localization-delocalization transition occurs at the condition 2 A — 3^ = 0, where ki = 0, i.e. the polymer 
delocalizes in the right half-plane. The localization length ^ ~ l/Zci becomes infinite at the transition. Notice that 
k2 (related to the chain's penetration length into the unfavorable solvent) stays finite at the transition. Using the 
relation for the energy per replica, ek = — with k2 given by Eq.(|l7|) we see that the energy at the localization- 
delocalization transition is finite and is equal to —2/3 A. This is due to the circumstance that in contrast to the 
localization of a QM particle in a potential well, where the energy at the transition is zero, the interaction of the 
stochastic copolymer with the interface is in fact nonlocal.. The solution for ^ < can be obtained by using Eq.(p^) 
with ek = —Dkf. 

Considering now the symmetric situation (Co = 0). From ( p^[]T7| ) we get 

k, = fVA, fe^^TA. (18) 
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Eq. ( p^ tells us that the ground state breaks down the symmetry of the Hamiltonian 7i„ with respect to the reflection 
X —X. An explanation can be given as follows: The condition Co = means that the fraction of the different 
types of monomers is on average balanced. However, the symmetry of the composition is true only for an ensemble of 
copolymers in the space of the quenched variable Ci^). The typical copolymer in such an ensemble contains an excess 
of species of one type, which causes the asymmetry of the localized state. It is remarkable that the Hamiltonian 7i„ 
contains information on this asymmetry and that our variational approach is able to reflect this as an asymmetry of 
the ground state. 

Let us now compare the extremal principle based on the functional J^{G) given by Eq.(^) with that of Rayleigh-Ritz. 
The latter is based on the minimization of the expectation value of the Hamiltonian iJ„, < k \ iJ„ \ k >, computed 
by using the trial function associated with (|l2). As in the case of localization of a particle in an asymmetric potential 
well considered above, the minimization of (1^) with respect to e, ki and ^2 completely coincides in the limit of 
large t with the Rayleigh-Ritz procedure. The difference between both methods appears to be due to the condition 
Ek = —Dk2. A restriction in the search of the bound state, which is imposed by this condition, results apparently in 
selecting the states breaking the reflection symmetry of the Hamiltonian Hn- This explains why in contrast to the 
Rayleigh-Ritz method our approach is successful in finding the localized state of the problem under consideration. 
While Eq. (|4|) still coincides with the stationarity condition of < A: | Hn \ k > with respect to ki , the left-hand side of 
Eq.(|l5|) differs from the derivative of < fc | iJ„ | k > with respect to fc2- The inspection of the condition dJ-/dk2 = 
shows that the latter provides the equality = < fc | Hn \ k > /n. This shows that the minimization machinery 
in both variational principle based on !F{G) and that of Rayleigh-Ritz is in general different, so that even in the 
stationary case the both extremal principles are not completely equivalent with each other. In the example of the 
localization of a QM particle in a potential well the both methods give, however, the same result. The equivalence of 
Eq. ( p^ with 9 < fc I Hn I fc > dki = guarantees that the energy is extremal with respect to fci . 

We have found that the energy associated with the asymmetric solution ( p^ ) is lower than the energy obtained 
from the minimization of Eq. (^3|) in the symmetric case fci = fc2 . The computation of the inverse Laplace transform 
of jl^ ) with respect to z gives the functional T at the extremum as: T{t) = — e'xp{Dnk2t) /t. Thus, the asymmetric 
solution (^8|) gives for all nonzero n lower values for J-'{t) than the symmetric one. 

To conclude, we have studied the adsorption of an ideal random copolymer chain at a selective interface. In terms of 
the replica method adsorption is related to the ground state of a replica Hamiltonian. To proceed we have introduced 
a novel variational principle for the Green's function. By using an appropriate trial Green's function as a possible 
candidate for the ground state we find the phase diagram of the localization-delocalization transition for a random 
copolymer chain on the interface. We predict that even in the case of a symmetric composition (Co = 0) the ground 
state is obtained to be asymmetric. This means that the ground state breaks down the reflection symmetry of the 
Hamiltonian. 
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